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Abstract 

We present preliminary results from a search for the flavor-changing neutral current decays B — ► 
JJtt+£- and B -> i<T* (892)^+^~ using a sample of 22.7 x 10 6 T(4S) -» SB decays collected with the 
BABAR detector at the PEP-II B Factory. We have reconstructed the following final states: B + — > 
K+^", B° -» K°£+^- (i^ -» vr+vr-), B+ -► K* + £ + £- (K*+ -► #°7T + ), and B° -» K*°l + £~ 
SZ • (K*° — » i^ + 7r _ ), where £ + £~ is either an e + e~ or n + n~ pair. We obtain the 90% C.L. upper limits 

23(B -» Ki + f) < 0.6 x 10" 6 and S(B -> iT7 + .T) < 2.5 x 10" 6 , close to the Standard Model 
predictions for these branching fractions. 
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1 Introduction 

The flavor-changing neutral current (FCNC) decays B — > K£ + £~ and B — > ET* (892)^ + £~, where 
£ is a charged lepton, are highly suppressed in the Standard Model, with branching fractions 
predicted to be of order 10 — 10 -6 . The dominant contributions arise at the one- loop level and 
involve afnf->s transition, so these decays are sensitive to the Cabibbo-Kobayashi-Maskawa 
factor V t * Vtb- The loop can involve either the emission and reabsorption of a virtual VF-boson, with 
the radiation of a virtual photon or Z (that subsequently materializes into the £ + £~ pair), or the 
emission of two virtual W bosons, producing the £ + £~ pair through a box diagram. Such processes 
are known as electroweak penguins. The simpler decays B — ► K* r y and B — ► X s j (where X s is any 
hadronic system with strangeness) have been observed jj], ||, providing the first evidence for the 
electromagnetic penguin amplitude. 

These rare decays are interesting not only as a probe of Standard Model loop effects, but 
also because their rates and kinematic distributions are sensitive to new, heavy particles that can 
appear virtually in the loop ||. Such heavy particles are predicted, for example, by supersymmetry 
(SUSY) models. 

Table 111 lists the predictions of a number of calculations based on the Standard Model. In 
such calculations, the short-distance physics that governs free-quark decay is incorporated into 
the Wilson coefficients Cf , Cg ff , and C\q in the Operator Product Expansion (n| of the effective 
Hamiltonian. This part of the calculation is relatively well understood in the Standard Model, 
except for nonperturbative contributions in the dilepton mass regions near the charmonium res- 
onances. There are, however, significant uncertainties in the predicted rates, which arise from 
strong (QCD) interactions among the quarks. Although these long-distance effects are difficult to 
calculate, they can be rigorously parametrized in terms of form factors, which are estimated with 
methods such as light-cone QCD sum rules or lattice QCD. From Table [l], we see that the rate 
for B — > K*£ + £~ is expected to be three to four times larger than that for B — » K£ + £~ . The 
amplitudes for the B — > K*£ + £~ modes include a pole at q 2 = "m 2 +e _ = 0, where the photon is on 
the mass shell. This pole gives a substantial contribution to the rate for kinematic configurations 
in which q 2 ~ and is particularly large in the B — > K*e + e~ mode. Measurements of B — > X^ 
constrain the magnitude of Cf , the coefficient of the electromagnetic penguin operator, but new 
physics could still affect the phase of Cf . 

We search for 5-meson decays in the following channels: B + — > K + £ + £~, B° — ► K°£ + £~ 
(K° -» 7T+7T-), B + -> K* + £ + £- (K*+ -> i^ s °7r+), B° -» K*°£ + £~ (K*° -» K + ir-) where £ + l' is 
either an e + e~ or a fi + fi~ pair. Throughout this paper, charge conjugate modes are implied. 
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Model K£ + £- K*e + e - K*fi + fi- 

LCSRQ 057^ 2^1, 1.9+^ 
LCSR |, | 1.4 

Quark Models §, 0.62 ±0.13 2.1 ±0.7 1.5 ± 0.6 
QCD Sum Rules |, § 0.3 1.0 

Table 1: Branching fraction predictions in various models within the Standard Model framework: 
light cone QCD sum rules (LCSR), quark models, and QCD sum rules. 



2 Detector Description and Data Samples 

The data used in the analysis were collected with the BABAR detector at the PEP-II storage ring 
at the Stanford Linear Accelerator Center. We analyzed the data taken in the 1999-2000 run, 
consisting of a 20.7 fb _ sample taken on the T(4S) resonance, as well as 2.6 fb _ taken at 
a center-of-mass energy 40 MeV below the T(4S) resonance peak to obtain a pure continuum 
sample. Continuum events include non-resonant e + e~ — > qq production, where q = u, d, s, or c. 
The on-resonance sample contains (22.7 ± 0.4) x 10 6 T(4S) — > BB events. 



The BABAR detector is described in detail elsewhere [11]. All components of the detector were 
used for this study. Of particular importance for this analysis are the charged-particle tracking 
system and the detectors used for particle identification. At radii between about 3 cm and 14 
cm, charged tracks are measured with high precision in a five-layer silicon vertex tracker (SVT). 
Tracking beyond the SVT is provided by the 40 layer drift chamber, which extends in radius from 
23.6 cm to 80.9 cm. Just outside the drift chamber is the DIRC, which is a Cherenkov ring-imaging 
particle identification system. Cherenkov light is produced by charged tracks as they pass through 
an array of 144 five- meter-long fused silica quartz bars. The Cherenkov light is transmitted to the 
instrumented end of the bars by total internal reflection, preserving the information on the angle 
of the light emission with respect to the track direction. The DIRC is used for kaon identification 
in this analysis and is essential to our background rejection. Electrons are identified using an 
electromagnetic calorimeter comprising 6580 thallium-doped Csl crystals. Muons are identified in 
the Instrumented Flux Return (IFR), in which resistive plate chambers (RPCs) are interleaved 
with the iron plates of the flux return. 

3 Overview of Analysis Method 

Our search for B — ► K^£ + £~ was guided by several key considerations. Given the predicted 
branching fractions and our reconstruction efficiencies (6% to 18%), we expect only a handful of 
observed signal events. Because the experimental signature is strong, however, it might be possible 
to establish a signal with a small number of events, as long as effective background suppression 
is achieved. But in working with small event samples, we must ensure that no bias is introduced 
through the event selection that could produce either an artificial enhancement of events in the sig- 
nal region, or an artificial suppression of events in the sidebands used to determine the background 
level. 

To avoid the possibility of bias, both the signal region and the sidebands used to measure the 
background levels are "blinded" during the process of defining the event selection criteria, so that 
statistical fluctuations in these regions cannot be exploited to artificially enhance a signal. Once 
the event selection is defined, we apply the cuts, unblind the data, and perform unbinned maximum 
likelihood fits to determine the signal and background yields. In these fits, the overall background 
level is determined from the data, with shapes derived from Monte Carlo but checked against 
control samples from the data. Although we define signal and sideband regions for the purpose of 
optimizing event selection cuts, the distinction between these regions is not relevant to the fit. 

To optimize the event selection without using either the signal region or the sidebands used 
for the final background determination, we have used (1) full GEANT 3 JL^] Monte Carlo for signal 
and backgrounds, (2) high statistics fast Monte Carlo for backgrounds, and (3) a continuum data 
sample collected at a center-of-mass energy below the T(4S). For simulation of signal events, we 
have used the event generator EvtGen [13] to implement a full form-factor parametrization of the 



decay amplitudes. The simulation provides the detailed kinematic correlations among the particles; 
in particular, the joint distribution in the Dalitz plot variables q 2 and lepton energy is simulated. 



Final-state photon radiation is simulated with the PHOTOS 14] package. The selection procedure 
was optimized using matrix elements from the model by Ali et al. |J, but the predictions from a 
number of other models were studied as well. About twenty thousand events were generated for 
each signal final state using a full GEANT 3 simulation. 

To validate the Monte Carlo modeling of both the signal and backgrounds, and to determine 
some of the systematic errors, we have used an extensive set of control samples from the data. 
The modeling of distributions characterizing signal events, and the associated cut efficiencies, were 
checked using B — > J/tpK and B — > J/vpK* decays (with J/ip — > £ + £~). These decays yield 
final states that are topologically identical to signal events (although they have a very narrow 
distribution in q 2 = m 2 +i _ ) . These modes are both serious backgrounds that must be suppressed 
and powerful control samples for understanding the signal efficiency. Another control sample is 
provided by reconstructed Kefj, and K*e^i combinations, which cannot result from B decays in 
the Standard Model. The Monte Carlo predictions for the background shapes associated with 
combinatoric backgrounds can be checked using these samples. Certain individual backgrounds, 
especially those that peak under or near the signal, are studied with GEANT 3 Monte Carlo samples 
of these decays. Finally, we defined a "large sideband" region in the data that is sufficiently far 
from the signal that we do not need it for determining the background level, but which allows us 
to check that expectations for background levels based on Monte Carlo are essentially confirmed in 
data prior to unblinding. 

The major sources of backgrounds are T(4S) —> BB and continuum events. More specifically, 
these include (1) B — » J/tpK^*' or B — > ip(2S)K^*' events; (2) generic BB backgrounds with 
either two real leptons or one real lepton and one hadron misidentified as a lepton (usually a 
muon); (3) background from continuum processes, especially cc events with a pair of D — > K^*'£ + v 
decays, or events with hadrons faking leptons; (4) very small but potentially serious contributions 
from a number of B decay modes with similar topology to the signal, such as B + — > D°ir + with 
D° —> K~ir + , in which hadrons misidentified as leptons can create a false peaking signal; and (5) 
BB or continuum events with photon conversions in the detector material. 

A large number of Monte Carlo samples were generated to study these backgrounds. We studied 
B — » J/t/jK^*' and B — > t/j(2S)K^*' events using about twenty thousand GEANT 3 events in each 
final state. Approximately 8.4 M GEANT 3 generic BB events and about 22 M continuum events 
were generated. In addition, we have developed a fast, parametrized Monte Carlo and generated 
220 M BB and 690 M continuum events, corresponding to 200 fb _1 of integrated luminosity. The 
realism of this simulation was enhanced by the fact that many key quantities are simulated using 
look-up tables with detailed information on charged-track efficiencies, particle ID efficiencies, and 
particle ID misidentification rates, all of which were measured in studies of data. 

4 Event Selection 

In the initial step of the event selection, we require that each event have at least four charged 
tracks and that the ratio of the second and zeroth Fox- Wolfram moments |15| be less than 0.5. 
This requirement provides a first suppression of continuum events, which have a more collimated 
("jet-like") event topology than BB events. We also require that the events contain two oppositely 
charged loosely identified leptons with laboratory-frame momenta p > 0.5 (1.0) GeV/c for electron 
(muon) candidates. Electrons are identified primarily on the basis of E/p, where E is the energy 

10 



measured in the Csl electromagnetic calorimeter. Muons are identified mainly by the number of 
nuclear interaction lengths penetrated by the charged track through the detector. 

The events are also required to lie within a large, rectangular region in the plane defined by two 
kinematic variables pi] ]: ttt-es > 5.0 GeV/c 2 and \AE\ < 0.8 GeV. Because signal events contain 
two B mesons and no additional particles, the energy of each B in the center-of-mass (CM) frame 
of the T(AS) must be equal to the e + or e~ beam energy in the CM frame. We define 



™ ES = K (£bc am ) 2 -(Ep:) 5 



AE = 2 \K + (P*) 2 - ^beam, 
i 

where E£ eam is the beam energy in the CM frame, p* is the CM momentum of particle i in the 
candidate 5-meson system, and mi is the mass of particle i. For signal events, the beam-energy- 
substituted B mass, tbeS) peaks at tub with a resolution of about 2.5 MeV. The quantity AE is 
used to determine whether a candidate system of particles has total energy consistent with the 
beam energy in the CM frame. 

We next apply a tighter set of particle identification requirements on both leptons and hadrons. 
The very tight selection criteria we apply to electrons give an efficiency of about 88% assuming well 
measured charged tracks, with a corresponding probability for a pion to fake an electron signature 
of about 0.15%. For our very tight cuts, the typical muon ID efficiency for momenta greater than 
about 1 GeV/c is 60%-70%, with a pion fake probability of about 2%. Electrons and positrons 
are required to pass a conversions veto, which suppresses e + e~ pairs that could have come from 
photon conversions in the detector material. For the B° — > K*°£ + £~ channels, we reconstruct the 
K*° in the K + ir~ final state, while the K* + in B + — > K* + £ + £~ is reconstructed in the K s tt + final 
state. The mass of the K* candidate is required to be within 75 MeV/c 2 of the mean K*(892) 
mass. The charged kaon daughters of the K* are required to be identified as kaons, and all charged 
tracks not identified as kaons are taken to be pion candidates. The kaon identification is based 
on the measurement of the Cherenkov angle from the DIRC for p > 0.6 GeV/c and on the energy- 
loss measurements (dE/dx) from the silicon vertex tracker and the drift chamber for momenta 
p < 0.6 GeV/c. The typical charged kaon ID efficiency is about 80% for well reconstructed tracks, 
with a pion fake probability of about 2%. K s candidates are reconstructed by combining two 
oppositely charged pions and requiring that the mass of the candidate be within 9.3 MeV/c 2 of 
the mean Kg mass, that the Kg candidate vertex probability be greater than 0.1%, and that the 
distance between the vertex and the candidate B vertex be greater than 0.1 cm. For each final 
state we select at most one combination of particles as a B signal candidate per event. If multiple 
candidates occur, we select the candidate with the greatest number of drift chamber and SVT hits 
on the charged tracks. 

Charmonium background, B — ► J/t/jK^-*' and ip(2S)K^*\ is suppressed by applying a veto in 
the AE vs. m£+£- plane. The boundaries of the vetos are shown in Fig. |l[ It is not sufficient to 
remove events with dilepton masses consistent with those of the J/ip or tp(2S), since bremsstrahlung 
or track mismeasurement can result in a large departure of the dilepton mass from the resonances. 
Charmonium events can pass this veto if one of the leptons and the kaon are misidentified as each 
other ("swapped"). There is also a fairly significant feed-up from the B — ► J/ipK and B — > ip{2S)K 
modes to the B — ► K*£ + £~~ channels. Both the swapped candidates and the feed-up contributions 
are explicitly vetoed. 
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Figure 1: The veto regions in the AE vs. rri£+£- plane that are populated by B — > J/ipK and 
B —> ^){2S)K events. The shaded areas are vetoed. (The same veto is applied to the K* modes.) 
For reference, the two horizontal lines bound the region in which most signal events are found. The 
charmonium vetos remove these backgrounds not only from the signal region, but also from the 
sideband region, simplifying the description of the background in the fits. 

As we have indicated, the charmonium decays B — > J/ipK^*' and B —>■ ip(2S)K^*' provide 
powerful control samples for checking the efficiency of our analysis cuts with events that have 
properties very similar to those of the signal. Figure |2| compares the AE distributions (absolutely 
normalized) of these charmonium samples in Monte Carlo vs. data. We find good agreement in 
both the normalization and the shape. 

Continuum background is suppressed using a four- variable Fisher discriminant, which is a linear 
combination of the input variables with optimized coefficients. The variables are the Fox- Wolfram 
moment; the cosine of the angle between the B candidate and the beam axis in the CM frame 
(cos#b); the cosine of the angle between the thrust axis of the candidate B meson daughter 
particles and that of the rest of the event in the CM frame; and the invariant mass tuki- The 
last variable helps discriminate against background from D semileptonic decays, for which m^e 
is below the D mass. BB combinatorial background is suppressed using a likelihood ratio, which 
combines candidate B and dilepton vertex probabilities; the significance of the dilepton separation 
in z (along the beam direction) ; cos 6b ', and the missing energy, E m iss , of the event in the CM 
frame. The variable E m [ ss provides the strongest discrimination against BB background: events 
with two semileptonic B decays often have more unobserved energy (due to neutrinos) than signal 
events. 

For the purpose of optimizing the event selection, a signal box is defined in the AE vs. w-es 
plane: 5.272 < m ES < 5.286 GeV/c 2 and -0.11 (-0.07) < AE < 0.05 GeV for the electron (muon) 
channels. Using normalizations from the GEANT 3 based Monte Carlo and shapes from the fast 
Monte Carlo, we perform extrapolations to estimate the number of expected signal events, S, and 
background events, B, in the signal region for a given set of event selection criteria. We select the 
criteria that minimize the probability that the backgrounds will fluctuate to a level at least as large 
as the predicted signal. The result of this optimization is very similar to an optimization based on 
maximizing S 2 /(S + B). 
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Figure 2: Comparison of AE shapes between data and Monte Carlo for the charmonium control 
samples, (a) B + -► J/VH e+ e -)K + , (b) B+ -► J/^(-> /J. + n~)K + , (c) 5° -» J/^(-» e+e")^* , 
(d) B° -> J/^H /z + aO**°, (e) 5° -+ J/^H e+e")^, (f) 5° -> J/^(-> m + M")^, (g) 
B+ -► J/^(-» e+e-)K*+(^ K§tt + ), and (h) 5+ -► J/^(-» m~V~)^* + H i^vr+). The points 
with error bars show the on-resonance data, and the solid histograms show the prediction of the 
charmonium Monte Carlo. All of the analysis selection criteria have been applied except for the 
charmonium veto, which is reversed. The normalization is absolute. 
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5 Unbinned Maximum Likelihood Fit 

To extract the signal and background yields, we perform an extended unbinned maximum likelihood 
fit in the AE vs. th-es plane. The fit region is defined by rriES > 5.2 GeV/c 2 and \AE\ < 0.25 GeV. 
There are three components in the fit: the signal, the BB background, and the continuum udsc 
background. The likelihood is given by 

g-( ra sig+«flfl+™cont) N 

C(n sig , n BB ,n cont ) = — ]^[ (n sig i=V sls + n BB P BB + n cont P 4 cont ), (1) 

where N is the total number of events in the fit region; i is an index over the events; and P slg , 
P , and p cont are the normalized probability density functions (which depend on ities and AE) 
for the signal, BB, and continuum background components, respectively. In the fits to the data, 
the shapes of these three components are fixed (as described below), so that the fit is used only to 
obtain the yields of these components. 

The signal shape, P slg , must describe a number of important properties of the signal events, in 
particular, the effects of photon radiation from the leptons. The distribution in ities is essentially 
a Gaussian centered at the P-meson mass, with a small radiative tail extending to lower values 
of tues- The width of this Gaussian, roughly 2.5 MeV/c 2 in both electron and muon channels, is 
determined primarily by the beam-energy spread. The AE distributions are much broader, and 
the effects of radiation are more pronounced. To describe these shapes and their correlations, we 



fit Monte Carlo signal events to a product of functions known as "Crystal Ball shapes" [16], one 
each for mES an d AE. Correlations are taken into account by allowing the parameters in the 
function describing m^s to depend on AE. The function describing AE is actually somewhat more 
complicated than the Crystal Ball shape in that a sum of two Gaussians is used, rather than a 
single Gaussian, to help describe the large radiative effects. 

Background events come from two main sources, T(4£) — » BB and continuum production of 
qq pairs. The BB background shape, P , is parametrized in AE as an exponential of a second 
order polynomial in AE multiplied by an Argus function JLTJ] in tties with the slope parameter 
taken as a second order polynomial in AE. The continuum shape, p cont ; is taken as the product 
of a first order polynomial in AE and an Argus function in 7tt.es- The parameters for the BB and 
continuum shapes are determined from fitting our 200 fb _1 sample of fast Monte Carlo events. The 
ability of the fast Monte Carlo to predict these background shapes is checked using the Kefj, and 
K*e\x control samples in the data (Figure |3|). In each final state, the shapes derived from Monte 
Carlo describe the data well. 

Backgrounds that peak under the signal in tties an d AE are not included in the fit, because we 
have suppressed their contributions to a level well under one event. 

In the fits we constrain the BB and continuum components to be greater than or equal to zero. 
For the signal yield we have a lower cut-off such that the total fit function is positive to avoid large 
negative yields in fits where there are no events consistent with signal. 

6 Results 

Before unblinding the signal region and the near sideband we estimated, based on Monte Carlo 
and off-resonance data, that there would be a total of 9.1 background events in the signal boxes 
in the eight modes we analyze. After unblinding, we used the sideband region to determine the 
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Figure 3: Fitted m-ES distributions in the (K, K*)e n^ control sample in data in slices of AE: 
AE < -0.11 GeV for (a) Ke ± ^ and (b) K*°e 



and (d) K*°e ± ^; 



± ^- 



±/i=F 



-0.11 < AE < 0.05 GeV for (c) Ke ± fi 
and AE > 0.05 GeV for (e) Ke ± ^ and (f) K*°e ± fi T . The dashed lines 
show the continuum component of the fit, and the solid lines show the sum of continuum and BB 
components. 

background in the signal region and found the total background to be 8.4 events. This is in good 
agreement with our Monte Carlo and off-resonance predictions. 

Figure |] shows the results of unblinding the signal and the near sideband regions. We emphasize 
that the determination of the signal yields is based on the fit, not on counting the number of events 
in the signal box. The fit results are shown in Fig. || and are summarized in Table |2|. By generating 
multiple toy Monte Carlo samples based on the fit probability density function (see Eq. [[]), we 
determine the signal yield such that 90% of the samples give a number of signal events larger than 
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Mode 


Signal 


90% C.L. 


Equiv. 


e 


(AB/B) e 


(AB/B) fit 


B/W-" 


B/10- e 






yield 


yield 


bkg. 


(%) 


(%) 


(%) 




90% C.L. 


B+ 


-» K+e+e- 


-0.2 


3.0 


0.6 


17.5 


±8.6 


±10.6 


-0.1 


0.9 


B+ 


-» K+pt+fT 


-0.2 


2.8 


0.4 


10.5 


±8.6 


±10.6 


-0.1 


1.3 


B°- 


-» K*°e+e- 


2.5 


6.7 


1.8 


10.2 


±10.5 


±10.6 


1.6 


5.0 
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-» K*VV~ 


-0.3 


3.6 
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±10.8 


±10.6 
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3.6 
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-► K°e+e- 
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4.7 
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-» K°n + ^~ 


0.0 
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9.6 
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±10.6 


0.0 
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B+ 


-» K*+e+e- 


0.1 
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±10.6 


0.1 


10.0 


B+ ■ 


-> K*+fi + fi- 
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±10.6 


3.3 
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0.7 


3.3 


B+ 


-» if* + eV T 


-0.4 


3.7 


0.7 


9.3 


±11.7 


±10.6 
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Table 2: Results from the fits to B — > K^*>t + i~ modes. The yield refers to the number of 
signal events given by the fit. The next column gives the 90% C.L. upper limit signal yield. 
The contribution of the background to the error on the signal yield, expressed as an equivalent 
background yield (see text), is reported in the next column, e is the signal efficiency, defined as the 
efficiency for the signal events to be in the fit region and does not include the branching fractions of 
secondary decays. The systematic error on the efficiency, (AB/B) e , and on the signal yield from the 
fit, (A£>/£>)fit! are listed. The last two columns give the central values of the branching fractions 
(B) and the preliminary upper limits on the individual modes at 90% C.L. 



that observed in the data. For the modes where we observed a negative yield, we place the 90% 
confidence level limit assuming zero events were observed. Table [2] lists an equivalent background 
yield determined as the square of the error on the signal yield in toy Monte Carlo simulations in 
which only the background components were included. Table || also includes the results for the 
lepton- number- violating decays B — > K^*'efX, where the signal efficiencies were determined from a 
phase space Monte Carlo simulation of these decays. 

There are two main categories of systematic uncertainties that affect the limits: uncertainties 
from the fit on the extracted number of signal events and uncertainties on the signal efficiency. 
Sources of systematic uncertainties are summarized in Table 0. 

The sources of systematic uncertainty affecting the extraction of the signal yield are the errors 
associated with the choice of the background and signal parametrization in the fit. To evaluate 
the systematic uncertainty on the background shape, we took the difference between the nominal 
result for the 90% C.L. event yield and the largest deviation found if one allows the background 
shape to be given either by a pure BB or a pure continuum shape. We assign an error of ±7% 
to the uncertainty in the knowledge of the background shape. The systematic error on the signal 
parametrization is established by varying the signal shape using the parameters obtained from 
fitting the charmonium control sample, and is of the order ±8%. 

In setting the upper limit, the systematic errors on the efficiency, (AB/B) e , and on the signal 
yield from the fit, (AB/B)fy, are added in quadrature, and the limit is increased by the correspond- 
ing factor. 
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Figure 4: Scatterplots of AE vs. toes after all analysis cuts in (a) B + — > if+e+e , (b) -B + — > 
K+n+n-, (c) 5° -► K*°e+e~, (d) 5° -» K*V + M", (e) #° -» ^e+e", (f) 5° -► Kg/x+/i~, (g) 
5+ -v K*+(^ Kgvr + )e + e", and (h) B+ -> #*+(-► i^vr + )/i + /i~- The small rectangle indicates 
the signal region, which is used only for optimizing event selection criteria. 
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Figure 5: Fitted ?ties distributions in the AE signal region for (a) B + — > K + e + e , (b) B + — * 
K+n + n~, (c) B° -» K*°e+e-, (d) 5° -► K*°n+n~, (e) 5° -► i^e+e", (f) 5° -» Kg^/i", 
(g) £+ -» #*+(-> i^vr+^+e", and (h) 5+ -> K* + (^ K^TT + )fi + fi- . The dashed lines show the 
continuum component of the fits, the dotted lines show the sum of continuum and BB components, 
and the solid lines show the sum of background and signal components of the fit. 
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Source of systematic uncertainty AB/B (%) 

Lepton tracking efficiency ±1.2 

Pion from K* tracking efficiency ±2.0 

Tracking efficiency for other tracks ±1.3 

Electron identification ±2.7 

Muon identification ±2.0 

Kaon and pion identification ±2.0 
Monte Carlo statistics for signal ±(3.0 to 5.0) 

Continuum suppression cut ±2.0 

BB suppression cut ±3.0 

K s selection ±4.0 

Modeling of signal kinematic distributions ±5.0 

Number of BB pairs ±1.6 

Background shapes ±7.0 

Signal shapes ±8.0 



Table 3: Sources of systematic uncertainties on the limits. 

7 Combined Results and Conclusions 

We average the electron and muon channels to determine the average branching fractions B(B — > 
K£ + £~) and B{B — > K*£ + £~). The modes are combined using the likelihood fit where a combined 
likelihood is formed as a product of the likelihoods of the individual modes. The fit extracts a 
combined signal event yield. For the averaging of modes with K* mesons the ratio of branching 
fractions B(B —> K*e + e~)/B(B — ► K*fi + /j~) from the model of Ali et al. (see Table Q) is used to 
weight the yield in the muon mode relative to the electron mode. The extracted yield corresponds 
to the electron mode. The combined fits give 

B(B -> K£ + £~ ) = (0.0 ± 0.3(stat.)) x 10~ 6 
B(B^K*£ + £-) = (0.7±l.l(stat.)) x 10~ 6 . 

As there is no evidence for a signal we evaluate the upper limits on these combined modes and 
obtain the preliminary results 



B{B- 


-» K£ + £-) < 0.6 x 10" 6 at 90% C.L. 


B(B- 


> K*£ + T) < 2.5 x 10~ 6 at 90% C.L. 



based on the analysis of a sample of 22.7 x 10 6 T(AS) -> BB decays in the BABAR 1999-2000 data 
set. These limits represent a significant improvement over previous results [|l^, 16, 2C] and are at 
the same level as the Standard Model predictions listed in Table [I]. 
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